Relativistic calculation of the two-photon decay rate of 
highly-excited ionic states 

Ulrich D. Jentschura^'^'G and Andrey Surzhykov^'^ 

^ Max- Planck- Institut fiir Kernphysik, 

66 ■ 

O ■ Postfach 103980, D- 69029 Heidelberg, Germany 

^ [ ^Institut fiir Theoretische Physik, Philosophenweg 16, 69120 Heidelberg 

^ Physikalisches Institut der Universitdt, 
^T) • Philosophenweg 12, 69120 Heidelberg 

(Dated: July 3, 2008) 

Abstract 

Based on quantum electrodynamics, we reexamine the two-photon decay of one-electron atoms. 
Special attention is paid to the calculation of the (two-photon) total decay rates which can be 
viewed as the imaginary part of the two-loop self-energy. We argue that our approach can easily 
be applied to the cases with a virtual state having an intermediate energy between the initial and 
the final state of the decay process leading, thus, to the resonance peaks in the two-photon energy 
distribution. In order to illustrate our approach, we obtain fully relativistic results, resolved into 
electric and magnetic multipole components, for the two-photon decay rates of the 35i/2 -^ l'S'1/2 
transition in neutral hydrogen as well as in various hydrogen-like ions. 
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I. INTRODUCTION 



Q 



Since the seminal work of Goppert-Mayer [Jj, two-photon decay rates of excited states in 



flflfla 



hydrogen-hke atoms and ions have been the subject of intense experimental |^, l3|, \4, Isj and 
theoretical [g, [Zl, |8|, l9|, [lO|, lul, ll^, ll3| studies. For many years, the investigations have dealt 
not only with the total decay rates but also with the energy and even angular distributions 
of the two emitted photons. By analyzing these (two-photon) properties, unique information 
has been obtained about the structural properties of one-electron systems including subtle 
relativistic effects as well as about the basic concepts of quantum physics such as, e.g., the 
entanglement. 

Even though large experimental and theoretical efforts have been undertaken in the past 
to understand various aspects of the two-photon decay of hydrogen-like atoms, the analysis 
of this process still raises a number of unresolved problems. One of these problems, which 
currently attracts a lot of interest, concerns those two-photon transition from highly excited 



states to the ground state which pass throu 
the initial state of the two-photon process 



A an intermediate state with a lower energy than 



10| and can alternatively decay to the ground- 



state via two (or more) sequential one-photon emissions. Such a transition leads to resonance 
peaks in the energy spectrum of the coherently emitted photons from the two-photon decays 
which are located at the energies corresponding to the (real) intermediate states. One of 
the most pronounced examples of such a situation is the 3Si/2 — > l'S'1/2 two-photon ElEl 
transition for which the differential (in energy) emission rate has singularities at energies 
corresponding to the 35*1/2 ~^ '^Pi/2 -^ l'S'1/2 and 35*1/2 ~^ 2P3/2 — > 15*1/2 cascades. A 
proper treatment of these singularities is obviously required for computing total decay rates 
obtained after an integration over the energies of the coherently emitted photons in the 
direct two-photon decay 35*1/2 -^ 15*1/2- 

During the last two decades, the theoretical treatment of the resonances in the energy 
distribution of the emitted photons has been discussed in a number of places. In general, 
the decay rate Fj of an initial state \i) in a hydrogenlike atom is the sum of a one-photon 
decay rate F^ and a two-photon contribution F 



I ) 



^(1) , r(2) 



F, = F^^ + Tr> . (1) 



Q 



The expression for F^^ as originally derived in [ij is easily seen to involve an integral over the 
energies of the emitted photons, the sum of which has to be equal to the energy difference of 



the initial and final states, and a summation over all possible intermediate, virtual states. In 
order to avoid problems with non-integrable singularities, the authors of [9] have attributed 

(2) 

Tl only to the so-called non-resonant intermediate transitions, in contrast to a summation 
over the complete intermediate-state spectrum. The non-resonant transitions are given 
by intermediate states of energy higher than the energy Ei of the initial state [9]] (the 
"resonant" intermediate states, which are involved in the one-photon cascade processes, are 
explicitly excluded from the sum over intermediate states). Based on this assumption, the 
non- resonant contribution for the 35*1/2 -^ l'S'1/2 two-photon transition was calculated as 
r35. = 8.2196 s~^ for hydrogen. Later, this result has been also confirmed in the calculations 



by Florescu and co-workers 10| who used a different method for the summation over the 
"non-resonant" states. 

Although the results presented in Refs. [9|, [lO(] are in mutual agreement, they are both 

(2) 
based on the interpretation of the two-photon decay rate Ti as a rate generated only by non- 
resonant intermediate levels. In our manuscript, we would like to propose an alternative way 
for the computation of the (two-photon) total decay rates which leads to a natural removal 
of the infinities otherwise introduced into the expression for the two-photon decay rate. 
We apply here a fully relativistic, quantum electrodynamical approach to re-investigate the 
two-photon decay of highly-excited states of hydrogen-like atoms, paying special attention 
to a careful handling of the resonances infinitesimally displaced from Feynman's photon 
integration contour (these singularities exactly correspond to the problematic "resonant" 
intermediate states). By making use of this approach, we obtain finite, physically sensible 
results for the decay rates of the two-photon 351/2 -^ 15i/2 transitions in neutral hydrogen 
as well as in the various hydrogen-like ions. Apart of the leading, electric dipole (ElEl) 
transition, we also discuss the contributions from the higher multipole components to the 
total decay rate. 

This paper is organized as follows: after a brief survey of the theoretical expressions used 
in our analysis (Sec. [TTl), we proceed by discussing the method of evaluation (Sec. IIIII) as well 
as the numerical results obtained for the differential as well as the total two-photon decay 
rates (Sec. HVT) . Conclusions are given in Sec. |Vl 



II. THEORY 



Within quantum electrodynamics, the (negative) imaginary part of the self-energy is just 



the r/2, where T is the decay width IJ, [l5|, [l6| . The one-photon decay T^^^ rate is obtained 



from the imaginary part of the one-loop self-energy, while the two-loop self-energy gives rise 
to the two-photon decay rate T^^'. Because the relativistic formulation of the two-loop self- 



energy problem has been discussed before in a number of places 17|, [iSj, we only mention 



here that by following a straightforward generalization of the standard procedure described 



for the non-relativistic framework in Refs. 



12 



I13I . lid ] , we obtain the following expression 



for the two-photon decay rate {h = c = tQ = 1), 



„2 



Fj = — limRe / dujiujiuj2 d^id^2 Si f{uji,uj 2) , (2) 

71 e^O J J 



where cui -|- u;2 = ciJmax = Ei — Ef with the initial and final state energies Ei and Ej, 
respectively. Sij is given by 

Ei - E^ - ui + ie 
^V Ei-Ep-iOi+ie 

{iJi\A2\ijp){i^p\A,\^f) \ 

Ei-Ep-uj2 + it )' ^ ' 

where in the second factor, the initial and the final state are exchanged, but the infinitesimal 
imaginary part in the denominators remains +ie (i.e., does not change sign). We here man- 
ifestly assume that ipiir) = i'niji^li{f) and ipjir) = 4'nfjf^if{f) are the well-known solutions 
of the Dirac Hamiltonian for a single electron in the standard representation, describing 
an electron bound to a point-like nucleus with charge number Z. For photons propagating 
with wave vector ki {i = 1, 2) and unit polarization vector n^. (fcj ■ Ux- = 0), moreover, the 
electron-photon interaction operator Aj in the transition amplitude (I3l) can be written in 
velocity gauge as: 

Ai = Aoa- ux^e'^^"- , (4) 

where Aq is a normalization factor, a are the standard Dirac matrices, and Aj = ±1 denotes 
the helicity, i.e. the spin projection of the photon onto the direction of propagation fcj. It is 



important to note that even though the electron-photon interaction operator (jl]) depends, 
of course, on the direction of the photon emission, one has to integrate over these directions 
in Eq. ([2]) in order to get the total decay rate. 

III. EVALUATION 

The summation over the intermediate states in the amplitude (l3l) runs over the complete 
one-particle spectrum \ip^) = \ipn„j„^i„), including a summation over the discrete part of 
the spectrum as well as an integration over the positive and negative-energy continuum 
of the Dirac spectrum. One has to use the full Dirac-Coulomb Green function — which 
is not known in closed analytic form — in order to perform this calculation consistently. In 
recent years, the Green's function method [19] has been widely applied for the analysis of the 



total two-photon decay rates as well as the photon-photon angular correlation functions |20 |. 
Various possibilities for the numerical implementation of the relativistic Green's function are 
known, among which we would like to mentionas (i) a well-known formulation in terms of 



Whittaker functions 



2l| and (ii) a Sturmian decomposition in terms of Laguerre polynomials 



as suggested by Hylton and Snyderman [22] 



We use the latter representation and apply the techniques of Racah's algebra to all 
spherical tensors and to the standard radial-angular representation of the wave functions, 
and to the Dirac-Coulomb Green function. For the interaction of electrons with the radiation 
field, the spherical tensor components are obtained from the known standard multipole 



expansion of the photon operator [see, e.g., Eq. (5) of Ref. |23|]. 



LMp 



1(P) 



where Af^j^ denotes the electric {p = 1) and magnetic {p = 0) multipole fields, respectively. 

IV. RESULTS 

The great advantage of the multipole decomposition ([5]) is that is allows us to study 
the contributions to the total (two-photon) decay rate from the various allowed multipole 
combinations. We use the integrand in the integral over uJi in Eq. ([3]) as a measure of 
the differential decay rate (where we can set explicitly e = for the differential rate). 
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FIG. 1: (Color online) Differential decay rate dT/dx for the dominant ElEl component of the 
3S — > 15 two-photon transitions in neutral hydrogen as well as in hydrogen-like ions with Z = 
1, 10, 20, 30, 35 and 40, where x is the energy sharing x = uji/{uji +uj2)- Relativistic wave functions 
are used for the initial, intermediate and final states, but the electron-photon interaction has been 
restricted to electric dipole term {ElEl term). The resonance peaks in the decay rate correspond 
to the resonant 35*1/2 ~^ '^Pi/2 ~^ l'5'i/2 ^-^id 35^/2 "*■ 2-P3/2 — > l'S'1/2 decay processes. 

The energy distributions of the two photons emitted in the 3S'i/2 — > 15'i/2 decay of neutral 
hydrogen and hydrogen-like ions are calculated as a function of the energy sharing parameter 
X = uji/{uji + 002). For an energy sharing in the range < x < 1, the contributions to the 
energy distribution from the ElEl and E1M2 multipole combinations are displayed in 
Figs. 1 and 2, respectively. As seen from these figures, the photon energy distributions for 
both multipole combinations exhibit sharp resonance peaks. As already mentioned, this 
behaviour is due to the fact that the summation in Eq. ([3]) includes also intermediate states 
|z/) having an energy E^, with Ei > E^ > Ej. However, the intermediate states contributing 
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FIG. 2: (Color online) Differential decay rate dT/dx for the E1M2 component of the 3S -^ IS two- 
photon transitions in neutral hydrogen as well as in hydrogen-like ions with Z = 1, 10, 20, 30, 35 and 
40 (x is the fraction of energy carried by the El photon). The E1M2 terin is treated relativistically. 
The resonances in the decay rate exclusively correspond to the 35i/2 -^ 2P3/2 -^ l'S'1/2 cascade, 
the 2Px/2 state does not contribute. 

to the peaks are not only defined by the (symmetry of) the initial \i) and final |/) states but 
also by the multipole components of the radiation field involved in the two-photon process 
and are different for ElEl as opposed to E1M2, and in addition, marked differences exist 
between the low-Z and the high-Z region. E.g., the fine-structure of the resonance in the 
ElEl energy spectrum grows with Z (the contributing states are 2Pi/2 and 2P3/2). By 
contrast, no splitting is observed — even for very heavy ions — for the ii^lM2 component of 
the 3S'i/2 -^ lSi/2 decay. Only one intermediate state, namely 2P3/2, is allowed for E1M2. 
Our treatment of the intermediate state resonance peaks in the integration over the 
photon energy is dictated by an accurate analysis of Eqs. ([2])^(I3]) as obtained form the 
relativistic two-loop self energy. The general structure of the contribution to Sif{u!i,u!2) 



generated by virtual states with intermediate energies Ei > E^ > Ef with resonance energy 

ujj. = Ei — Ey or ujr = Ey — Ef is as follows: 

Sif{uJi,UJ2) ^— — + - ^ — —2- (6) 

UOr-^ + lt {uJr-UJ + le) 

The integration of the first term can be carried out using the Dirac prescription (see, for 
example, Ref. [l^). The second term of Eq. ([6]) can be treated using the formula 

1 . o . 



limRe/dcj =^ -, (7) 

where we used an appropriate scaling of the photon energy integration variable in order to 
map the integration region to the interval (0, 1). It is important to note that Eq. ([7]) holds 
strictly for < cu^ < 1, but the limit is not approached uniformly |12l . [l3|; i.e., it would be 
forbidden to exchange the sequence of the limit e ^ with the integration over uj. As usual 
in quantum electrodynamic processes, all regulators have to be kept up until the very end 
of the calculation. 

With these preparations, it is easy now to integrate over the photon energies (see Table 
1 for the 35*1/2 ~* l'S'1/2 process). As seen from the table, the cross sections for the ElEl, 
E1M2 and MlMl components components scale with the nuclear charge as Z^, Z^^ and 
Z^^, respectively. As expected, this scaling behaviour is similar to the Z-scaling of the 
multipole components in the 251/2 — ^ l'S'1/2 transition. 

Furthermore, as seen from Table 1 and as implied by the non-uniform convergence of 
the integrals, the intermediate states with the energies E^ lying between the energies of the 
initial and the final states give a finite contribution to the two-photon decay rate. For the 
electric dipole (ElEl) transition in a neutral hydrogen atom, e.g., a proper treatment of 
the intermediate 2Pi/2 and 2P3/2 states leads to the decay rate of FgJ = 2.08 s~^ which is 
in agreement with the result of nonrelativistic calculations reported in Ref. [12j]. However, 
when comparing our prediction with the theoretical data by Cresser and co-workers [9|] a 
large discrepancy by about a factor of 4 is observed. The occurrence of the discrepancy 
is natural because the problematic virtual states with intermediate energies are treated 

CI 

differently in [9|. 
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2.08 


1.19 


6.13 
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2.03 


1.18 


6.13 


10 


1.98 


1.16 


6.14 


15 


1.94 


1.14 


6.16 


20 


1.90 


1.12 


6.20 
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1.84 


1.08 


6.24 


30 


1.79 


1.03 


6.30 


35 


1.67 


0.96 


6.39 


40 


1.60 


0.86 


6.50 



TABLE I: Contributions from different combinations of multipoles to the integrated decay rate 
r(2)^ in units of radians per second. Relativistic calculations liave been performed for different 
hydrogen-like ions. 

V. CONCLUSIONS 

The two-photon decay of hydrogen- like ions has been re-investigated within the framework 
of relativistic quantum electrodynamics. Starting from first principles of this theory, we 
treat the total (two-photon) decay rate as the imaginary part of the relativistic two-loop 
self-energy. The great advantage of this approach, which has its roots in field theory, is that 
it provides a simple and efficient route to handle the potentially problematic cases of those 
two-photon transitions from an excited into the ground state which pass intermediate states 
that can otherwise also be reached in one-photon cascades from the initial to the final states. 
We found that those states with energies that lie between the energy of the initial and the 
final states, contribute a finite correction to the total two-photon decay rate. Taking into 
account this correction, we calculate the rates for the 3S'i/2 — > 15'i/2 two-photon decay of 
neutral hydrogen as well as hydrogen-like ions. Our results are in a good agreement with 
nonrelativistic calculations for low Z (see Refs. |12l . [l3|) but show a significant deviation 



from the data by Cresser and co-workers 



9| 



Our quantum electrodynamics approach, as discussed in the present paper, opens a way 
for a systematic theoretical analysis of the simultaneous, coherent two-photon emission from 
one-electron (and many-electron) atomic systems, even in cases where problematic interme- 
diate states with an energy between the initial and final states give rise to resonance peaks 
in the photon energy distributions. We stress here that a conceivable alternative approach 
to the removal of the formal infinities generated by the intermediate "cascade" states, which 
is based on the explicit removal of these states from the sums over v and p in Eq. ([2]), gives 
rise to a number of conceptual problems, including gauge-noninvariance with respect to 
length and velocity gauges [131]. Our approach is manifestly gauge invariant and also avoids 
problems connected with the identification of the infinitesimal parts ie in the propagator 
denominators in Eq. ([2]) with partial or total decay rates of the intermediate states: the e 
parameters are free parameters which approach zero after all other operations, including the 
integrations over the photon energies, have been performed. This operation leads to a finite 



result and corresponds, as explained in Ref. 13|, to a partial removal of the problematic 



intermediate states from the sum over all virtual states involved in the two-photon process, 
albeit in a fully gauge-invariant manner. 

In addition to its relevance for atomic physics, our approach may have a significant im- 
pact for astrophysical studies where a detailed knowledge of the (properties of) two-photon 
transitions is highly required for the analysis of cosmological hydrogen and helium recom- 
bination. The contribution of two-photon processes to the recombination history represents 



an issue which has recently attracted substantial theoretical interest |24| . 
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